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Abstract 

Source-sink systems are metapopulations of patches that can be of variable habitat 
' quality. They can be seen as graphs, where vertices represent the patches, and the weighted 

oriented edges give the probability of dispersal from one patch to another. We consider 
either finite or source-transitive graphs, i.e., graphs that are identical when viewed from 
a(ny) source. We assume stochastic, individual-based, density-independent reproduction 
and dispersal. 

By studying the path of a single random disperser, we are able to display simple criteria 
for persistence, either necessary and sufficient, or just sufficient. In case of persistence, 
we characterize the growth rate of the population as well as the asymptotic occupancy 
frequencies of the line of ascent of a random survivor. Our method allows to decouple the 
roles of reproduction and dispersal. Finally, we extend our results to the case of periodic or 
random environments, where some habitats can have variable growth rates, autocorrclated 
in space and possibly in time. 

In the whole manuscript, special attention is given to the example of regular graphs 
where each pair of adjacent sources is separated by the same number of identical sinks. 
In the case of a periodic and random environment, we also display examples where all 
patches are sinks when forbidding dispersal but the metapopulation survives with positive 
probability in the presence of dispersal, as previously known for a two-patch mean-field 
model with parent- independent dispersal |15) . 
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1 Introduction 

Model. We consider a stochastic, individual-based model of spatially structured population 
dynamics. The spatial structure is a metapopulation of patches that can be of different habitat 
qualities. We label hy i = 1, . . . , K the patches so that the model can be described by a labeled 
finite graph with weighted oriented edges. Vertices represent the patches, an oriented edge 
from vertex i to vertex j bears a weight dij equal to the probability of dispersal from patch i 
to patch j. 

We assume a simple asexual life cycle with discrete non-overlapping generations and no 
density-dependence. At each generation, as a result of survival and reproduction, all indi- 
viduals, independently from one another, leave to the next generation a random number of 
individuals, called offspring, whose mean number depends on the habitat quality. Immedi- 
ately after local growth, each individual from the new generation migrates independently, from 
patch i to patch j with probability dij. No mortality is assumed during dispersal, since it can 
be encompassed in the growth phase. Reproduction, survival and dispersal probabilities are 
assumed not to depend on local densities. 
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The mean per capita number of offspring in patch i wih be denoted by rrii, and can be seen as 
a proxy for habitat quahty. If > 1, the population may persist without dispersal and we say 
that patch i is a source. Conversely, if mj < 1, the population would die out without dispersal 
and we say that patch i is a sink. It will be convenient to assume that mi > m2 > • • • > 
and more relevant to consider the case when mi > 1 > mx- Even in the presence of sinks, the 
metapopulation might persist thanks to local growth on sources replenishing sinks by dispersal. 
In the second part of the paper, we will also assume that environment is variable and denote by 
mi{w) the mean offspring in a patch of type i when the environment is in state w. Finally, we 
extend our results to source-sink metapopulations on transitive graphs. We also want to make 
the observation that our approach applies to more general multitype branching processes, in 
the sense that dispersal behaviours of siblings might not be independent. 

Two natural examples. Let us describe two examples with two possible habitat qualities, 
one source type and one sink type, that will be treated throughout the paper: 



First, we are interested in the simple case with two patches, patch 1 with mean offspring M, 
and patch 2 with mean offspring m. In this case, we will always use the simplified notation 
p = di2 and q = d2i (see Figure [I]). 



Figure 1: Two patches of different qualities. The filled circle is a source and the empty circle 
is a sink. The arrow labels are the dispersal probabilities. 



Second, we will consider the case when each source is only connected to sinks and two adjacent 
sources are separated by a line of n identical sinks. An example of such graph is the cyclic 
finite graph with one source and n sinks, or two sources connected by n sinks, or an infinite 
array with period n (see Figure E])... 



Outline of the paper. A random disperser is a walker on the graph which follows the 
dispersal stochastic scheme. We display here a very simple persistence criterion involving the 
mean growth rates in each type of habitat and the times spent in each habitat, between two 
visits of a source, by a single random disperser. We specialize this result for the examples 
given above. We also give the asymptotic fraction of time spent in each patch by an individual 
taken at random in the surviving population. We find that this distribution never equals the 
stationary distribution of the single disperser, except in the case when all habitat qualities are 
identical. 



M := mi > 1 



m := m2 = ■ ■ ■ = mx < 1. 




p 
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Figure 2: Three examples where sources are only connected to sinks and two adjacent sources 
are connected by n sinks; a) cyclic graph with one source, n = 7; b) two sources, n = 3; c) 
infinite periodic array, n = 2, arrow labels (not represented) are also assumed periodic. 

Addressing those questions is much more difficult when the environment is variable, i.e., 
when each habitat type can have a different mean growth rate depending on the state of the 
environment. Applying our techniques yields partial results, viz. when the environment is 
periodic or given by an ergodic sequence of random variables. In these two cases, we prove 
that the population can survive in sinks only, in the sense that without dispersing, it would die 
out in each patch, a result which was previously known |15j in the case of large populations on 
an island model (parent-independent dispersal). 

Finally, we extend naturally our approach to a wide class of infinite graphs that are called 
transitive, but with a finite set of habitat qualities. As phrased in [21], a graph is transitive 
if it globally "looks the same" from any vertex (including labels and edges). The vertices of 
the graph are now labeled by a countable set D and dpq gives the dispersal probability from 
patch P G -D to patch Q € D. The type of patch P, denoted by i = j{P) G {1, • • • K}, gives 
its quality and the mean number of offspring per indivual in patch P is equal to m^i^py In this 
paper, we will assume that the graph is only source-transitive, that is, transition probabilities 
and habitat qualities seen from any patch of type 1 are identical. We specify the mathematical 
definition in Section [5l In particular, a practical example is given by sources with the same 
quality connected by corridors of identical sinks and of the same length (see Figure [6] for an 
example). A transitive graph could also be an (infinite) chessboard where whites are sinks 
and blacks are sources, the square lattice 1? where sources have coordinates of type (n, n) 
(diagonal) or of type (ra, 0) (horizontal array), and so on (see Figures [3] and S]) .. . 



Background. In ecology and conservation biology, an important issue is to understand the 
effects of spatial heterogeneity on population persistence, in particular in the context of global 
climate change and wildlife preservation management. In particular, recent papers have focused 
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a) 



b) . 




Figure 3: Two examples of source-transitive graphs with a finite motif ; a) the chessboard; b) 
a square grid where four-degree vertices are sources separated by n sinks (here n = 1). 

on the factors (dispersal, temporal autocorrelation) that can allow for persistence even in sink 
metapopulations, see e.g. [151 13 [IBJ and references therein. 

In mathematics, the stochastic process describing the size of a density- independent, individual- 
based, stochastic population is called a branching process [1, 3j. Here, we consider multitype 
branching processes, where the type of an individual is given by its dwelling patch. When 
the environment varies, the population can be described by multitype branching processes in 
varying or random environment, see [21 \W[ [20] for mathematical references. 

(i) 

The number of individuals located in patch i in generation n is denoted by Zn ■ The 
process Z = {Zi'\ i = !,••• n > 0) is a multitype Galton-Watson process. It is known 
from the mathematical literature [U [3] that eithei0 the population becomes extinct or it grows 
exponentially. More specifically, we see that ruidij is equal to the mean number of offspring of 
an individual living in patch i which will land into patch j in one time step, and therefore we 
call mean offspring matrix the matrix A defined as 



The maximal eigenvalue (see e.g. [19]) of A is called Malthusian growth rate, or simply growth 
rate. Indeed [31 [l], if /) < 1, the metapopulation dies out with probability 1, and if p > 1, the 
metapopulation can survive with positive probability, in which case 



where W is multidimensional finite and positiv( 



^we exclude in this paper the degenerate case when the size of population is constant with probability 1 
■^under some integrability assumptions on the offspring numbers and the strong irreducibility (or the primi- 
tiveness) of the matrix A. These conditions will be fulfilled here. 



A := {ruidij : i,j = 1,--- ,K). 




Z 



P' 
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Figure 4: Two examples of source-transitive graphs embedded in 7? with an infinite motif ; a) 
a diagonal of sources; b) a horizontal array of sources. 



Method. We call random dispcrser a single individual who moves on the graph at discrete 
time steps following the dispersal probabilities. In other words, if X„ denotes the position of 
such a random disperser after n time steps, then is the Markov chain with transition 

probabilities 

F{Xn+l=j\Xn = i)=dij. 

The goal of this paper is to display persistence criteria, along with results regarding the asymp- 
totic growth rate and the asymptotic fraction of time spent in each patch (by an individual 
taken at random in the surviving population) in terms of some characteristics of the motion 
of the random dispcrser. In contrast with the method involving the maximal eigenvalue of 
the mean offspring matrix, this one can yield quite simple, interpretable and partially explicit 
criteria. In addition, these criteria decouple the effect of dispersal and reproduction on popu- 
lation survival. This approach is still valid when the graph is an infinite transitive graph with 
finite motif, and partially valid even with infinite motifs. In a number of remarks, we will also 
provide sufficient conditions for survival which are explicit. 



Notation 


Interpretation 


vertex i 


patch 


oriented edge 
with weight dij 


probability of dispersal from patch i to patch j 


rui 


mean number of offspring in a patch with habitat type i 


M = mi 


mean growth rate in a source habitat (case mi > 1 > m2 = ■ ■ ■ = mx) 


m = 1712 


mean growth rate in a sink habitat (case mi > 1 > m2 = • • • = mx) 


P = di2 


probability of dispersal from the source to the sink (case of 2 patches) 


q = d2i 


probability of dispersal from the sink to the source (case of 2 patches) 


X 


random walk on the graph following the dispersal probabilities 



Table 1: Notation. 
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We assume that at least one loop-edge has a positive weight, that is, there is at least one 
patch in which the probability of staying put is nonzero. We further make the assumption 
that any patch can be visited with positive probability by the random disperser whatever its 
initial position. The previous two assumptions make the Markov chain (X„;n > 0) irreducible 
and aperiodic and the matrix D (and A) primitive (or strongly irreducible). Actually, our 
results could be extended to this framework and even to the general case by lumping states 
together into classes but we will not develop this point here. We do not consider here the 
degenerate case when the per capita offspring number is constant and a.s. equal to 1 in patch 
1. Finally, we assume that the offspring number Ni of any individual living in patch i satisfies 
E(iVj log"*" A'^j) < oo (finite A log A moment). This guarantees the convergence of Zn/p^ to a 
non degenerate r.v. W which is positive on the survival event. 

2 A first result on global persistence 
2.1 General case 

We now give a criterion for metapopulation persistence in terms of the random disperser X. 
For that purpose, we assume from now on that the random disperser starts in patch 1 [Xq = 1) 
and we denote by T the first return time of the random disperser into patch 1, 

T := min{n > 1 : A„ = 1}. 



Theorem 2.1 The population persists with positive probability iff 

miE(nfj-/mxJ > 1. 

We observe that 

where Sn{i) '■= #{1 <k<T—l: X^ = i} is the time spent in habitat i by time T, that is, the 
number of times the random disperser has visited patches of habitat type i until time T. 

Proof. Let a (ancestor) be some individual placed at time in patch 1. Define Yi as the 
number of offspring of a staying put in patch 1. Now for any integer n > 2, let Yn denote the 
number of descendants of a at generation n living in patch 1 and whose ancestors at generations 
1, 2, . . . , n — 1 have all lived outside patch 1. Then set 

n>l 

that can be seen as the total number of descendants of a who live in patch 1 for the first time 
in their lineage (except a). 

In the theory of random trees, the set of individuals belonging to one of the Yn individuals 
for some n, is called a stopping line (see e.g. [5j). By the branching property of our tree, the 
total numbers of descendants of each of the individuals of this stopping line who live in patch 1 
for the first time in their lineage, are independent and all follow the same law as Y. Therefore, 
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the total number of descendants of a living in patch 1 is finite iff the branching process with 
offspring number distributed as Y is finite, which is equivalent to the a.s. extinction of Y. 
The bottomline is that the population persists with positive probability iff K{Y) > 1. Indeed, 
we have excluded the case when the per capita offspring is a.s. constant equal to one in patch 1. 

Let us compute K{Y). We first note that for every i = 1, . . . ,K, 

K 

E{z(})) = ^E{zi^l,)mjdji, 

3=1 

(i) 

where Zn denotes the number of individuals located in patch i at generation n. We prove 

(i) 

easily by induction that the number of individuals Yn in patch i at generation n which have 
avoided patch 1 at generations k = 1,2, - ■ ■ ,n — 1 satisfies 

K 

J=2 ji,...,jn-ie{2,...K} 

As Yn = Y^^\ we get 

ji,...,jr^-ie{2,...K} 
= miE{\T=nmximx2 ■ ■ ■ ?Tix„-2'Tix„_i)- 

Adding that Y = J2n>i ^^^^ ' i^dive 

E{Y) = mi ^ E(lr=n"7,xi • ■■mx^_i) = miE(mxi • ■■mxr-i)- 

n>l 

This yields the result. □ 
2.2 Case of two habitat types 

Let us focus now on the special case when there are 2 habitat types and the source is solely 
connected to sinks : 

M := mi > 1, m := m2 = ■ ■ ■ = uik < 1- 

We denote by 

K 

p = Y 

the probability of dispersing for an individual living in patch 1. The per capita mean offspring 
number sent out from a source at each generation is Mp. Let a be the time of first visit of a 
sink by the random disperser 

a := inf{n > : X„ 7^ 1}, 
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so that (7 is a geometric random variable with success probabihty p. Next, let S denote the 
waiting time (after a) before the random disperser visits a source (this source might or might 
not be the initial source patch Xq) 

S := inf{n > : X^+„ = 1}. 

The duration S can be seen as the time spent in sinks between two consecutive visits of sources. 
In the two-vertex case, and in this case only, it is a geometric random variable with success 
probability q. The previous theorem reads as follows. 

Proposition 2.2 The population persists with positive probability iff 

M{l-p) + eMp>l, (1) 
where e is the depleting rate due to the sink habitat in the graph, defined as 

e := E{m^) = ^m¥(5 = k). 

k>l 

To see that, we use the first transition of the random disperser to get E,(IlJ^^mXi) — (1 ~P) + 
pE{m^). 

Remark 1 // the average time spent in sinks has 

M - 1 



E{S) < 



Mp{l — m) ' 



then the population persists with positive probability. Indeed, the mapping / : x i— )• f{x) = E{x^) 
is convex so 

e = f{m) > 1 + /'(l)(m - 1) = 1 - (1 - m)E{S) > 1 - = ^ ~ ~ 

which yields eMp + M(l — p) > 1. 
Remark 2 Criterion (JT]) is equivalent to 

EiWm^) > 1. (2) 
Both criteria are equivalent formulations of Theorem \2.1\ 



In contrast with the computation of a maximal eigenvalue, our criterion allows to decouple 
the roles of the mean offspring numbers M and m from that of the metapopulation structure 
itself (encompassed in S and p, or S and a), that is, of reproduction and survival. 

One can check in some simple cases that criterion ([T]) or the condition that the maximal 
eigenvalue p oi A exceeds unity are equivalent. For example, in the case when there are only 
one source and one sink [K = 2) and two vertices. Here the mean offspring matrix A is 




M{l-p) Mp \ 
mq m{l — q) J ' 
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The characteristic polynomial C of this square matrix is 

C{x) = (M(l — p) — x){m{l — q) — x) — Mmpq. 

Either M(l — p) > 1 and the population living in the source ensures the persistence. Or 
M{1 — p) < 1 and the quadratic polynomial is convex and has non negative derivative at 1. 
Thus, its leading eigenvalue is greater than 1 iff C(l) < 0, which reads 

Mp 1 — m(l — q) 

1 — M{1 — p) mq 

We recover ([T]) since here S is geometric with success probability q, which yields 

e = yq{l- qf-^jn^ = -. 

k>l ^ ^' 

Notice that even in this simple case where ^4 is a 2 x 2 matrix, the computation of the lead- 
ing eigenvalue is cumbersome, and we have used a trick to explicitly specify the persistence 
criterion. 



2.3 Example with pipes of identical sinks 

Assume that the source is a vertex of degree 2 in the graph, connected to a left sink and a 
right sink. The probability of staying put on a source is still 1 — p, the probability of dis- 
persing onto a left sink is pL, and the probability of dispersing onto a right sink is pR (so 
that L + R = 1). The sinks form a pipeline of n adjacent sinks linking adjacent sources. 
The probability of staying put on a sink is always s, the probability of dispersing from a 
sink onto one of its two neighboring sinks is r in the left-to-right direction of the pipe, and / 
in the right-to- left direction of the pipe (so that q = l + r = l — s). See Figure [5] for an example. 




Figure 5: A pipeline where n identical sinks (n = 7) connect the source to itself. 



This example will be directly extended in the last Section to infinite transitive graphs, 
where pipelines of n sinks periodically connect sources (see Figure ED . 
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pL 




1-p 




p 




Figure 6: A periodic pipeline where adjacent sources are separated by n identical sinks (n = 3). 

We can compute exactly the depleting rate e of the above kind of source-transitive graph. 
Let X > 1 > iJ, he the two ordered solutions to 



mrx^ — — ms)x + mZ = a; > 0. 



Then 



and 



A + /X = 



1 — ms 



mr 



Proposition 2.3 The depleting rate e is equal to 

Remark 3 In the 'chessboard' case (n = l), we recover 

(1 — s)m 
1 — ras 

In the case when dispersal is isotropic {I = r), we get 

_ A" - /i" + A - /X 

Notice that in both previous cases, the depleting rate does not depend on L or R. In the case 
of one single source and a large number of sinks (n oo), we get 

e = LX~^ + Rfi. 

In the case of one single source and isotropic displacement, we then get e = X~^ = /x. 

Proof. Consider a random walk y on {0, 1, . . . , n + 1}, with displacement at each timestep 
being —1 with probability I, with probability s, and +1 with probability r. Let Tj denote 
the first hitting time of i by Y, and set T := min(ro, r„+i) as well as 

Ofc :=E(m^ \Yo = k). 

Returning to the random disperser X on the graph, it is easily seen that 

E{m^ I is a left neighbour) = ai, 

while 

K{m^ I Xfj is a right neighbour) = a„, 
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so that 



e = Lai + Ra. 



Computations of ai and a„ rely on the fohowing recurrence relationship 



k G {l,...,n} 



with boundary conditions ag = On+i = 1- This relation is obtained easily by considering the 



3 Growth rate and habitat occupation frequencies 
3.1 General case 

If u denotes an individual in generation n, we define i/fc(u) as the patch occupied by the 
ancestor of u in generation k < n. Then, for every i G {1, . . . , K}, 



is the occupancy frequency of patch i by the ancestral line of u. We further denote by an 
individual chosen randomly in the surviving population at generation n. We will see that the 
dispersal history of U„, as described by (Fj(U„) : i = 1, . . . ,K) can be very different from 
that of a random disperser, since the ancestors of surviving individuals have better chance of 
having spent more time in sources than in sinks. 

It is known that the growth rate p of the metapopulation is equal to the maximal eigenvalue 
of A (see [HEIII?]). Moreover, the asymptotic occupancy frequencies of U„ are deterministic 
and can be expressed as the product of the right and left maximal eigenvectors associated to p. 
We refer to |13 1 I14 ^ [TT]. to Theorem 3.1 and 3.2 in f9] in continuous time and to [17j in discrete 
time. In this section, we want to give an alternative characterization to these quantities in 
terms of the random disperser and show an application. 

For that purpose, we use the transition matrix D of the random disperser X on the graph. 
We denote by the occupancy frequency of patch i by the random disperser X by time n 



By assumptioio, the random disperser has a a stochastic equilibrium on V, that we denote 
hy u = {ui : i = 1, .. . ,K), which is the unique positive solution to uD = u. By the ergodic 
theorem, we also know that with probability 1, 



A typical single disperser will therefore occupy patch i with asymptotic frequency Ui. This 
may not be the case of the ancestors of surviving individuals, whose paths must have favoured 
source patches. There is a trade-off between the pay-off in terms of fitness, gained by visiting 
source patches, and the cost in terms of likelihood, paid by deviating from the typical dispersal 
behaviour. This trade-off is particularly obvious if we consider the case of a perfectly sterile 

''the random walk X on the graph is irreducible and aperiodic 



first transition of the walk X. 



□ 




n 



FiiXn) ■.= -#{0<k<n:Xk=i}. 




12 



patch where the mean offspring is zero. In this case, the path followed by the ancestors of 
a surviving individual will necessarily have avoided this patch. Nevertheless, the asymptotic 
occupancy of this patch by a random disperser must be nonzero by assumption. 

There is a way of quantifying both the cost and pay-off of deviating from the typical dispersal 
behaviour, that is, of having asymptotic occupancy frequencies / = (/j : z = 1, . . . , K), where 
/ is a given element of the set J- of non-negative frequencies on the graph 



f = {fi:i = l,...,K):fi>0, J2 = ! 



i=l,...,K 



First, the probability that a random disperser has occupancy frequencies close to some given 
/ by time n decreases exponentially with n with rate /(/), which can thus be interpreted as 
the cost of the /-occupancy scheme: 

/(/):= sup I Yl fi^og{vi/{vDy^:v»0^, (3) 

where v » denotes a positive vector, that is, Vi > for each i = 1, . . . , K. Indeed, large 
deviations theory [6j \7\ ensures that for any e ^ 1, as n — )■ oo, we have 

nfi - e < F^iXn) < /, + e for alH = 1, • • • , X exp(-n/(/)). (4) 

We refer the reader to forthcoming Section [A] for a more rigorous formulation. Taking v = u 
in ([3]) shows that /(/) is of course always non-negative. This function is also convex. When 
f = u, one can easily check that each partial derivative of / is zero and it can be proved that 
the supremum in ([3]) is attained for v = u, so that I{u) = 0. This was indeed expected, since 
f = u is the natural occupancy scheme of the random disperser. 

Second, the reproductive pay-off of / can be defined as the fitness of a non-random disperser 
with given /-occupancy scheme, that is 



i=l,...,K 

Indeed, the total size of a population of individuals all adopting this dispersal behaviour can 
be seen to grow like 

Ug.mf^ = exp(ni?(/)). 

Thus the cost (in terms of likelihood) for a population to follow some occupancy scheme is 
quantified by / and the reproductive pay off by R. The best strategy is to have an asymptotic 
occupancy frequency ip which maximizes the difference R — I. If this optimum is positive, then 
the population survives with positive probability. In addition, the ancestral line of a randomly 
chosen surviving individual will have visited patch i with frequency ipi. These results are stated 
below. The last assertion indicates that this optimal occupancy scheme 93 = = 1 . . . ,K) 
is always different from the natural occupancy scheme u of one single random disperser, except 
when all habitat types have the same quality. 
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Theorem 3.1 The growth rate p of the metapopulation is given by 

log(p) = max{i?(/)-/(/) i/eJ-}. 

In addition, for any patch i = 1,. . . ,K, conditional on the population being alive at time n, 
the occupancy frequency of patch i by the ancestral line of a randomly chosen individual U„ in 
the surviving population at time n, converges to ipi in probability 



Fi(U„) — > (pi, 
where the frequency vector ip & is uniquely characterized by 

log{p) = Ri^)-Ii^). 

The occupancy frequency 93 coincides with the stationary distribution u of X (if and) only if 

mi = m2 = . . . = mK- 

In the same vein, we refer to Theorem 3.3 in [S] for a description of the lineage of surviving 
individuals for multitype branching processes in continuous time. 

The proof is deferred to Section El Strong irreducibility guaranteed by the assumptions 
given in the Introduction is useful (only) for the last part of the Theorem. 

3.2 Example with two patches 

Let us apply this result for the simple case of two patches. Recall that p is the probability to 
go from patch 1 to patch 2 and q the probability to go from patch 2 to patch 1 (see Figured]). 
We compute /(/i, f2) which amounts to finding the maximum of 

u G [0, 1] ^ /i log(n/(n(l - p - q) + q)) + f2 log((l - u)/{u{p + q - I) + I - q)). 

We search the roots of its first derivative : 

\u (\—p — q)u + qJ \ 1 — u [p + q — l)u + \ — q J 

Now we focus on the case p + q = 1, assumption which is made in [15] and makes computations 
easier (the general case requires to solve a third degree polynomial equation). Roughly speaking, 
this corresponds to letting a fraction q of the population live in patch 1 (and p = ^ — q live in 
patch 2). Then we check directly from the last displayed equation that 

u = fi, l-u = /2 = l-/i 

and 

/(/l,/2) = /llog(/l/(Z) + /2log(/2/(l -g)). 

To compute max {ii(/) — /(/) : / € J-}, we just need to maximize 

/i G [0, 1] ^ /i log(mi) + (1 - /i) log(m2) - /i log(/i/(?) - (1 - /i) log((l - /i)/(l - q)). 
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We search again the roots of the first derivative 



log(mi/m2) - log(/i/g) - 1 + log((l - /i)/(l -q)) + l = 0, 



which yields the habitat occupation frequency of patch 1 : 



1-ipi 




(1 - q)m2 



qrrii 



We get 



= 1-(P2 



qrrii + (1 - q)m2 



and also recover the expected result p = qrrii + (1 — q)m2. 

4 Fluctuating environments 

4.1 General setting 

We now enrich our model with a fluctuating environment. The environment is embodied by a 
certain value w which belongs to a finite set of states. We assume that the environment affects 
simultaneously all patches, but not necessarily in the same way. We keep on assuming a simple 
asexual life cycle with discrete non-overlapping generations and no density-dependence. Now 
the environment is assumed to affect the reproduction scheme, but not the dispersal behaviour. 
Specifically, at each time step, conditional on the state w of the environment, individuals 
reproduce independently according to some distribution which depends on the habitat type of 
their dwelling patch. We denote by mi{w) the mean offspring number of individuals dwelling 
in patch i when the environment is in state w. 

Except in the last subsection, we will assume that the environment alternates periodically 
at each time step between two states. Actually, the same method would allow to deal with 
any finite number of environmental states varying periodically, which allows for the modeling 
of seasonal effects. In contrast with [15] who deal with large populations sending out at each 
generation given fractions of their offspring to each patch (parent-independent migration, also 
known as island model), we consider discrete populations with stochastically varying size, 
individually stochastic dispersal behaviours, and allow for metapopulation structure of the 
stepping-stone type. In this more general setting, we come to the same conclusion that in 
fiuctuating environments, the population may survive in sinks only. This means that in the 
absence of dispersal, the population in each patch would become extinct. But positive dispersal 
probabilities may allow for global survival. 

We call ei and 62 the two possible states of the environment, so now we have 2K habitat 
qualities mi{ej), for i = 1, . . . , K and j = 1, 2. Now the process Z = {Zn\i = 1, . . . , K;n > 
0) is no longer a homogeneous Markov chain but a multitype branching process in varying 
environment. However, restricting the observation of the metapopulation to times when the 
environment is in the same state allows to adapt the arguments of the previous section. Indeed, 
(^2^^ i = 1, . . . , i^T; n > 0) is a multitype branching process with mean offspring matrix A with 
generic element 




i,j = 1,... 



K. 



k=l,...,K 
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In the following subsection, we treat the case of two patches and determine the global persis- 
tence criterion. We then handle the general case using the random disperser X. 

4.2 Example with two patches and two periodic environments 

For simplicity, even if the environment is now variable, the two patches are still called respec- 
tively the source (patch 1) and the sink (patch 2). The mean number of offspring in the source 
are denoted by Mi = mi(ei) and M2 = mi(e2). In the sink, they are denoted by mi = m2(ei) 
and m2 = 1712(62) ■ 

Theorem 4.1 A necessary and sufficient condition for global persistence is 

MiM2(l - pf + (Mim2 + miM2)pq + mim2(l - qf > min (2, 1 + MiM2mim2(l -p- qf) ■ 

Remark 4 It is easy to find examples where both patches are sinks when forbidding dispersal 

but the metapopulation survives with positive probability in the presence of dispersal. Indeed, 
in the absence of dispersal, each patch is a sink if (and only if) M1M2 < 1 and mim2 < 1. 
Assuming for example that p = q = 1/2 and mi = m2 = m, the global survival criterion 
becomes M1M2 + m(Mi + M2) + > 4, which holds as soon as m{Mi + M2) > 4. 

Proof. The mean offspring matrix of (^2n; > 0) is given by 

a(l,l) = MiM2{l - pf + Mim2pq 

a(l,2) = Mim2p{l - q) + MiM2{l - p)p 

a(2,l) = miM2q{l - p) + mim2{l - q)q 

a(2,2) = mim2{l — q)^ + miM2qp. 

The maximum eigenvector of the matrix A = {a{i,j) : 1 < < 2) is the largest root of the 
polynomial 

x^ - (a(l, 1) + a(2, 2))x + a(l, l)a(2, 2) - a(l, 2)a(2, 1). 
So it is less than 1 iff 

a(l, 1) + a(2, 2) + ^{a{l, 1) - a(2, 2)f + 4a(l, 2)a(2, 1) < 2 

Then the criterion for a.s. extinction of {Z2n : n G N) is 

a(l, 1) + a(2, 2) < 2 and (a(l, 1) - a(2, 2))^ + 4a(l, 2)a(2, 1) < (2 - a(l, 1) - a(2, 2)f. 

The second inequality becomes a(l, 1) + a(2, 2) < 1 + a(l, l)a(2, 2) — a(l, 2)a(2, 1), which gives 

MiM2{l — p)^ + Mim2pq + mim2(l — q)'^ + miM2qp < 1 + MiM2mim2{l — p — q)^ ■ 

This completes the proof. □ 
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4.3 Global persistence for more than two patches 

Here, we extend the previous result to the case of a general, finite graph. We want to state 
a global survival criterion which generalizes Theorem 12.11 to periodic environments. Assume 
again that the random disperser starts at time in patch 1 and set T the first even time when 
the random disperser goes back to habitat 1 

T := min{n > 1 : X„ = 1 and n is even}. 

Theorem 4.2 The population persists with positive probability iff 

mlE(^fJ-/mx,(u;^)) > 1, 

where the sequence {wi : i > 1) can take one of the two values (ei, 62, ei, . . .) or (62,61,62,...), 
depending on the initial environment. 

This theorem can be proved by a direct adaptation of the proof of Theorem 12.11 replacing Zn 
with Z2n- 



4.4 Rate of growth and habitat occupation frequency 

The generalization to periodic environments of the results of the previous section can be 
achieved by changing the state-space {1, . . . ,K} of the random disperser to the state-space 
of oriented edges of the graph, i.e., ordered pairs of vertices 

£:={!,. ..,Ky. 

Denote by B the transition matrix of the Markov chain {X2n, X2n+i',n > 0), which indeed 
takes values in £. Then denote by T the set of frequencies indexed by £ 

T:={{fE,Ee£): fE>0, YfE = l}, 

and define the new cost function J : — )• M as 

/(/) := sup I fElog{vE/{vB)E) :v»o\, 

where v denotes a non-negative vector indexed by £, such that v » 0, that is, ve > for all 
E G £. Also define the new pay-off function : — ?> M as 

R(.f) ■= Y fElog{mi{ei)mj{e2)). 
E={i,j)e£ 

We can also provide an expression of / in terms of the entropy function using Theorem 3.1.13 
in |6]. The generalization of Theorem 13.11 can be stated as follows. 

Theorem 4.3 The growth rate p of the metapopulation is given by 

21og(p) = max W) -/(/):/ G J-} 
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In addition, for any patch i G {1, . . . ,K}, conditional on the population being alive at time n, 
the frequencies of occupation of patch i by the ancestral line of a randomly chosen individual 
JJn in the surviving population at time n, converges in probability : 

iG{l,...,K} 

where the vector {fij : {i,j) E £) is characterized by 

logip) = Ri^)-Ii^). 
The proof follows that of Theorem 13.11 with now 

E(|Z2„+il) = E (n[Lomx,,mx,,+0 = E (ug,Uf^,rn,{ei f-"'(^m,{e2f"'^^^^ 
where \Zn\ is the total number of individuals in source patches at generation n and 
S^^Hi) = #{k<n: X^k = i}, S^^^i) = #{k < n : X^k+i = i). 

Example with two patches. We extend the computations of the previous section to periodic 
environments. We focus on the case p + q = 1, which ensures that the mean fraction of 
individuals in patch 1 is equal to q. We recall that this is an assumption made in [15j . 
Let us first compute /(/) := svly>{Y1iE&s fE^og{vE/{vB)E) ■ ^ ^} '^here 

B(i,j){k,i) = djkdki = CfeC/, i,j,k,l = 1,2, 

where ei = q = 1 — p, €2 = p = 1 — q- Then, assuming that XIegS ~ 

vBk,i = ^ VEeuei = ek€i, 



so that 



and 



/(/) = sup I ^ fE log (^^) ■.V»0,Y,VE = l\ 

1^ VE.eEj 

Theni?(/)-/(/) = Y^j^^^ fElog{mEi{ei)mE2{e2)eEieEj fE)- Writing /(2, 2) = l-T,Ee£-{2,2) fE 
and using that the partial derivatives of i? — / equal in 99, we have : 

log {■mE^{ei)mE2{e2)^Ei^E2) + 1 - log(<^£;) - [log (m2(ei)m2(e2)p^)) + 1 - log(v9(2, 2))] = 
which gives the habitat occupation frequencies 



'PE = (^EieE2'mEi{ei)mE2{e2) 



eEieE2'mEiiei)mE2{e2) 



EdS 
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We can also now deduce the growth rate p. Using that 



eEif^E2mEi{ei)mE2ie2) = q^MiM2 + qpMim2 + pqmiM2 + p^mim2, Je = I, 

we get 

R{(p) - I{(f) = - \og{q^MiM2 + qpMim2 + pqmiM2 + p'^mim2) 

and 

p= {qMi+pmif/'^{qM2 + pm2Y''^ . 
This is the same growth rate as the one computed in |15] . 



4.5 Some comments on random environments 

A more natural way of modeling fluctuating environment in ecology is to assume random 
rather than periodic environment. The approach developed for periodic environments cannot 
be extended to random environments directly. Indeed, since the environment affects the whole 
metapopulation simultaneously, the randomness of environments correlates reproduction suc- 
cess in different patches. The process {Zn\i = 1, . . . , K;n > 0) counting the population sizes 
on each patch is now a multitype branching process in random environment (MBPRE). 

Let us denote by A(w) the mean offspring matrix (involving dispersal) in environment w. 
Specifically, the generic element aij{w) of A{w) is the mean offspring number of a typical 
individual dwelling in patch i sent out to patch j by dispersal, when the environment is w, so 
that 

Aij{w) = mi{w) dij. 

We will now assume that the state-space of environments is finite and that the sequence 
{wn : n > 0) of environment states through time is a stationary, ergodic sequence, possibly 
autocorrelated, in the sense that the states need not be independent. Under this assumption, 
it is proved @ in [5] that the limit 7 of the sequence 

- log II A{Wn)A{Wn-l) . . . A{wo) \\ 

n 

exists with probability 1 and is deterministic, where || . || denotes the maximum row sum of the 
matrix. This is interesting to us because it is further showrl in O [TBI [20] , that the extinction 
criterion and the growth rate of this MBPRE are respectively given by the sign and the value 
of 7, more specifically, 7 = log p. 

Unfortunately, this does not give a very explicit condition for global persistence. But again 
using the random disperser, we can give some sufficient conditions for survival. For simplicity, 
we turn our attention to the example of two patches and two environments ei and 62. At any 
time step, the probability that the environment is in state ei is denoted by z/ € (0, 1) (so the 
probability that the environment is in state 62 is l — u). We show again that the population may 
survive in sinks only. As in the case of periodic environments, the mean number of offspring in 
the first patch (the source) is denoted by Mi = mi(ei) and M2 = m2{ei). In the second patch 
(sink), they are denoted by mi = m2(ei) and m2 = ""1-2 (62). 

''under the further assumption E(log"'" || A{wo) \\) < 00, where expectation is taken w.r.t. the environment 
^again under some further assumptions, see Annex. 
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For the sake of simplicity, wc state the results for two special cases of ergodic environments: 
the case when the states are independent at each time step (with probability v and 1 — i/), and 
the case when the sequence is a Markov chain. We denote by a. the transition from e\ to 62 
and by /3 the transition from 62 to e\ . Then it is well-known that v = (3/{a + f3). The case of 
independent environments is recovered when a + P = 1. Note that as soon as a + ^ 7^ 1, the 
sequence of environment states is auto-correlated. 

Proposition 4.4 We have the following lower bound for the growth rate of the metapopulation. 

log(/9) > 1^ log(Mi) + log(m2) + va logipq) + - a) log(l -p) -|- (1 - - /3) log(l - q). 

Remark 5 Observe that this lower bound does not depend on M2 and mi. Again one can 
display examples where both patches are sinks when forbidding dispersal but the metapopulation 
survives with positive probability in the presence of dispersal. In the absence of dispersal, each 
patch is a sink if (and only if) M^M2~'^ < 1 and m'{m\~" < 1. Actually one can manage 
to keep 7 > while M^"*^ < 1, mi < 1 and m2 < 1, for example with M2 small and Mi 
large for some fixed p, q, mi, m2. This corresponds to 62 being a catastrophic environment in 
the source patch but the population survives in patch 2 when a catastrophe occurs. 

Proof. We consider only the subpopulation avoiding patch 1 when the environment is equal 

to 62- This means that this population reproduces with mean offspring number Mi(l — p) in 
patch 1 while the environment is ei. Each time the environment 62 occurs, we consider the 
part of this population which has dispersed to patch 2. This corresponds to a mean offspring 
number of Mip. This population then stays in patch 2 and reproduces with mean offspring 
number m2(l — q) until the environment is again equal to ei. We then consider the part of this 
population which goes back to patch 1. This corresponds to a mean offspring number of m2q- 
Thus the patch of the ancestors of the individuals we keep is equal to 1 (resp. 2) if it lived 
in environment ei (resp. 62). Then at time n, the mean size of the population we consider is 
equal to 



where Ni{n) (i € {1, 2}) is the number of times before generation n when the environment is 
equal to and Nij{n) {i,j € {1,2}) is the number of one-step transitions of the environment 
from Cj to type ej until time n. By ergodicity, we know that these quantities have deterministic 
frequencies asymptotically. In the case of a Markovian sequence of environments, as n — >■ 00, 



Nii{n) ^ v{\ — a)n, Ni2{n) ~ van, N2i{n) ~ (1 — i')fin, N22{n) = (1 — z/)(l — P)n. 

Using the growth of this particular part of the whole population directly gives us a lower bound 
for 7: 

7 > z/logMi-|-(l— i^) logm2H-i^alogp+(l— logq+u{l—a) log(l— p)-|-(l— :^)(1— ^) log(l— g). 




Ni{n) ~ vn, 



N2{n) ~ (1 - u)n, 



and 



Noticing that (1 — i')^ = vol completes the proof. 



□ 
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Remark 6 We could improve these results by considering more sophisticated strategies. For 
example, we could consider the suhpopulation which stays in patch 1 if (and only if) the number 
of consecutive catastrophes is less than k and then optimize over k. 

Actually, the proof relies on a stochastic coupling. Roughly speaking, the subpopulation we 
consider avoids the bad patches at the bad times and follows a (one type) branching process in 
random environment 611,612,621 and 622 respectively with stationary probabilities z^(l — a), va, 
(1 — u)l3 and (1 — — /3)and mean offspring Mi(l — p), Mip, m2q and 7712(1 — q). 
Observe also that we can derive a lower bound using the permanent of the mean matrix M of 
the MBPRE from Proposition 2 in But this lower bound is not relevant for understanding 
the survival event in sinks only. 

5 Metapopulation on infinite transitive graphs 

We now turn our attention to infinite graphs with a special property called transitivity. We 
label the patch of the graph with a countable set D. Each patch P ^ D has a type i = j{P) E N 
which gives its quality. That is, the mean number of offspring in patch P is equal to nij^py 
The transitivity assumption means here that the graph, including dispersal probabilities and 
habitat types, looks the same from any patch of type 1. 

Let us provide the reader with some examples. An infinite linear periodic array of patches 
is a transitive graph (see Figure ?? for example). This means that two consecutive sources 
are separated by the same sequence of habitat qualities. The patches can then be labeled by 
i E D = Z and the type of patch i is equal to f{i), for some integers N,K > and a function 
/ : Z — >■ [1,K] such that f{i + N) = f{i) for every i S Z. The chessboard is a transitive graph 
(see Figure [7]) . Star sources with 2d pipelines of sinks form a transitive graph (see Figure [8] for 



Let us now specify mathematically these definitions. We consider here graphs whose vertices 
are patches which can be labeled by P € D, the oriented edges from P to Q are weighed by 
dpQ and the type of vertex P is denoted by j{P). A mapping T of the graph is called an 
isomorphism if it conserves the types of the vertices as well as the weights of the oriented 
edges: T is a bijection of N such that for all P, Q € N, 



The associated equivalence relation ~ between the patches of the graph is defined by 

P ~ P' iff there exists an isomorphism T of the graph such that r(P) = P'. 

The class of a patch P is defined as Cl{P) = {P' : P' ^ P}. Every patch of this class has the 
type of P. The graph is then called source-transitive if all patches of type 1 belong to a single 
class. The collection of the distinct classes {Cl{i) : i ^ V) form a partition of the patches of 
the graph. Such subsets V of patches are called motifs. With a slight abuse of notation, the 
transition probabilities on a motif V are denoted by (dpg : P ^V,Q &V) where 



d = 2). 



dT{P)T{Q) — dpQ 



3{T{P))=3{P). 





Q'eCliQ) 



which does not depend on the element P E Cl{P). 
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Figure 7: a) The chessboard (with periodic arrow labels - not represented) is a source-transitive 
graph that can be collapsed into: b) a two- vertex graph (loop edges are not represented). 



The initial graph can be seen as a family of copies of a motif properly connected. Observe that 
not all transitive graphs have a finite motif. For example, the cases illustrated by Figure [3] can 
be collapsed into a finite motif, but those given in Figure H] (a diagonal or a horizontal array 
of sources in Z^) cannot. Indeed, Figure [7] shows that the infinite chessboard in can be 
collapsed into a 2-vertex graph and Figure [8] shows that star sources with 2d pipelines can be 
collapsed into a finite motif with one source and d pipelines of sinks. We can easily extend the 



Notation 


Interpretation 


vertex P 


patch 


set V of vertices 


motif 


type of P j{P) = i 


habitat quality of patch P is i 


oriented edge with weight dpQ 


probability of dispersal from patch P to patch Q 


rrii 


mean number of offspring in a patch with habitat type i 



Table 2: Notation. 



results of the previous section to infinite transitive graphs by considering the random disperser 
Xn on the graph which follows the dispersal probabilities [dpg : P,Q ^ D"^). In that purpose, 
we assume that the random disperser starts in a patch of type 1 {j{Xq) = 1) and we denote by 
T the first return time of the random disperser into a (possibly different) patch of type 1. 

T := mm{n > 1 : j(X„) = 1}. 

Theorem 5.1 The population persists with positive probability iff 

miE(nfj/m,(^^)) > 1. 
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Figure 8: a) Star sources with 2d pipelines of sinks (and periodic arrow labels - not represented) 
form a source-transitive graph (here d = 2) that can be collapsed into: b) a single source with 
d pipelines (loop edges are not represented). 

The number of individuals located in patch P in generation n is denoted by Nn ■ Collapsing 
the graph into some motif V , we denote by 

P'eCi{P) 

the total number of individuals in some habitat of the class of patch P £ V. In the case when 

(p) 

the motif is finite, Z = [Zk ' , P eV;n>0) is a multitype Galton-Watson process. We see 
that mj(^p-jdpQ is equal to the mean number of offspring of an individual living in patch P 
which will land into patch Q in one time step, and therefore we call mean offspring matrix 

A := {m^(^P^dpQ : P,Q£V), 

Then, everything happens as if the metapopulation evolves on a finite graph corresponding to 
the motif. We then derive the following counterpart of Theorem 13.11 We denote now by D the 
transition matrix of the Markov chain j(X„), and still use for / = (/p : P 

I{f) := sup{ E fplog{vp/{vD)p) : t; » 0}, 

and 

pgv 
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Theorem 5.2 We assume that the graph is transitive and has a finite motif. 
The growth rate p of the metapopulation is given by 



log(p) = max I R{f) - I{f) : ^ /p = 1, /p > I . 

In addition, for any P , conditional on the population being alive at time n, the occupancy 
frequency of a patch with type j{P) by the ancestral line of a randomly chosen individual U„ 
in the surviving population at time n, converges to in probability 



Fi{lJn) — > <fi, 

where the vector belongs to T and is characterized by 

log{p) = Riip)-I{ip). 

The occupancy frequency ip coincides with the stationary distribution u of X only i/mj(p) = 
"T'j{P') for all patches P,P'. 

The results obtained for the periodic and random environments can be derived similarly to 
the infinite transitive graphs. 



A Proof of Theorem 13.11 

Recall that is an irreducible, aperiodic Markov chain with transition matrix D. As the 

state space {1, . . . ,K} is finite, is strongly irreducible and there exists no € N such that 

-D"° >> (all positive entries). 

Now the real number /(/) is the cost for the habitat occupation frequencies associated with 
the random walk X to equal /. It gives the geometric decrease of the probability that the 
portion of time spent in habitat i until generation n is close to fi : 

. . . , /x) = lim hm -- logP(/fe - e < Ffc(X„) < A + e). (5) 

This result holds for any irreducible aperiodic random walk X when / is finite. Indeed, thanks 
to Sanov's theorem (see e.g. [H Theorem 3.1.6 page 62], / is convex and continuous. The 
function / is called the rate function associated to the path of the random walk X. 

Finally, the fact that the offspring distribution Ni for an individual living in patch i satisfies 
^{Ni log^ Ni) < oo ensures that for all i,j = l,...,K, 

E{Z( log+ Z[^'' I = 1> Zll"^ = for k ^ i) < oo. 

Then Zn/p^ converges to a non degenerate variable, see e.g. [3l Chapter 5, section 6, Theorem 
1]- 
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Expression of the growth rate p and habitat occupation frequencies. We start with 



one individual in patch 1. Recalling that for every i = 1, . . . ,K , E(Z^*^) = Y^^=i ^(^n-i)^j'^j 
we get by induction 



K 



K 



i0 = l, 3n=i, 
l<il,-Jn-l<A- 

This yields 

E(|z„|) = E(n^roWO 

Denoting by the number of visits of the random disperser X in patch i = 1, . . .V, : 

Sn{i) = #{k<n-l: 3{Xk) = i] = nFi{Xn), 

we deduce 

E(|Z„|) = E(^nf=im^ 

1. 

= / exp(n V fi log(mi))P(Fi(X„) G d/i, . . . , G d/^). 

Using ([5]) and Laplace method, we get 

log(E(|Z„|)i/") max{5^ /, log(m,) - /(/i, ...,/;,):/ e J^}. 

This proves the first part of the result. 

The maximum oi h := R — I is reached for a unique frequency ip, which means that there 
is a unique such that log(p) = R{^p) — This is the object of Proposition lA.ll Moreover 

the partial derivatives of /i at 93 are zero. As fx = ^ — fi — ■ ■ ■ fn-i, for every 1 < i < K — 1, 

d 

log(mi) - log(mA') - -^I\f=^ = 0. 

If if = p, then the partial derivatives of / at 99 are zero. This can be directly computed or 
deduced from / > and I{p) = 0. This ensures that for every 1 < i < K — 1, log(mj) — 
log{mK) = 0, i.e. mi = m2 = . . . = mx- This proves the third part of the theorem. 

Finally, let us prove that the habitat occupation frequencies of a typical individual are given 
by the vector if. Let e > and i G {1, . . . , K}. The individuals alive in generation n are labeled 
by Ufc, A; = 1, . . . , Z„. 

Following the first part of the proof, 



k=l 



K 

/ l|/,-^,|>eexp(n V/aog(mi))P(Fi(X„) G dfu . . .,FK{Xn) G d/x). 
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Using again ^ and the Laplace method, we get 



-logE( J^l|^^(,^)_^^l>,) Ci.^ 



n 

k=l 



with Ci,e = maK{J2f=i fi log(mi) - . . . , /j^) : / G J", \fi-ipi\> e}. The uniqueness of the 
argmax if ensures that Ci^e < Q.o- Moreover the growth rate Cj^o is equal to log(p) and 



1 1 

-logE(^l|^^(,^)_^^l>,) - -logp" Ci,, - a,o < 0. 

k=l 

Then 



fc = l 

In other words Ylk=i '^\Fi{ui,)-ip,\>e/ goes to in probability. Adding that Zn ~ Wp"^ a.s. as 
n — > oo and {VF > 0} = {Vn G N, > 0} a.s. ensures that 

Zn 

'z? 



ra— ^oo 

A:=l 



in probability. By dominated convergence, 



/ ^ Zn 

\ " fc=l 



fc)-Vi|>f 



n— Q 



Then, conditionally on Z„ > 0, denoting by Un an individual chosen uniformly in generation 

P(|Fi(U„) -V9i| > e, Z„>0)"^0. 
This proves that i<i(U„) in probability and completes the proof. 

Study of / and uniqueness of argmax R — I. The supremum I defined by 

K 

I{f) = /(/i, . . . , Ik) := sup{ /, \og{uj/{uD)j) : n G M^, » 0}, 

i=i 

is reached for a unique unit positive vector. This means that there exists a unique u{f) = 
{ui, . . . , uk) such that 

K 

Hf) = Yl ^Og{Uj{f)/{u{f)D),), ui{f) + ---+ UK{f) = 1, > 0, . . . , ukU) > 0. 

Indeed this vector u{f) realizes a maximum for u G M^,u >> and thus satisfies for j = 
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This equation characterizes u, see Exercise IV. 9 page 46 in [7], which ensures that u{f) is 
uniquely defined. Note also that if fi is the stationary distribution, uj = fj satisfies this equa- 
tion since {fD)i = /j, so that / = 0. 



Proposition A.l There exists a unique (p ^ J- such that \og{p) = R{ip) — 

Proof. We observe that / i-^ u{f) can be extended from T to [0, oo)^ \ {0} and satisfy 
on [0, oo)^ \ {0} by setting 



K 

u{f) = u{f/ \\ f \\), where |[/||=^/j. 



i=l 

Then R — I can be extended to [0, oo)^ \ {0} with 

{R - I){f) = ^^i^ log{mMf)D)j/u,{f)), 

j=l Z^k=l Jk 

and {R - I) (A/) = {R - /)(/) for every A G (0, oo). 

Consider a vector / which realizes the maximum of R — I and does not belong to the 
boundary of [0, oo]^. Then the partial derivatives are zero and for every i = I, . . . , K, 



1 



log{mi{uD)i/ui) + ^ fj 



dfi ^""-^^i _ dfi '"j 

{uD)j Uj 



Y.f=i fj log{mj{uD)j/uj 



■13 



[Ef=i fk? 



Using dU]) we get 



. wS'^^^j _ fj , duk _ ^ duk fjdkj _ duk fk 

j^^ {uD), - j^^ {uD), ^ dU ~ ^ dU j^, (uD), - A. dn Uk ■ 

so that 



, , / m / ^ T.f=ifj'^og{mj{uD)j/uj) 
log[mi[uD)i/Ui) = — 



Y,k=i fk 

Observe that the right hand side does not depend on i, so that for every i = 1, . . . , K , 

{uD')i = aui, 

where Dj- = djirrii and a is a positive constant. Then u is left eigenvector of D' with positive 
entries. Moreover D' is strongly irreducible since D is strongly irreducible and mj > for 
every i by assumption. Now Perron-Probenius theory ensures that there is a unique left positive 



D is strongly irreducible if there exists no > 1 such that all the coefhcients of D"" are positive. This property 
is a consequence of irreducibility and aperiodicity. 
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eigenvector u of D' such that Yli^i Ui = 1 (see e.g. [12] 0. Moreover, fohowing the hterature 
on large deviations (see e.g. [7j), ([6|) reads 

i=l 

with Dj- = Ujdji/{uD)i. Here again D" is strongly irreducible since D is strongly irreducible 
and both u and uD are positive vectors. Using again Perron-Frobenius theory guarantees the 
uniqueness of the solution / such that fi = ^- This ensures the uniqueness of the argmax 

of i? — I in the interior of [0, oo)^. We complete the proof by adding that there is at least one 
argmax in the interior of [0, oo)''^ since we recall that the frequency occupation is the product 
of the right and left eigenvectors of A associate to p, which are both positive. If (fi and ip2 
realize the max of R — I, the concavity of this function ensures that the segment [ipi, (p2], which 
is contradiction with the uniqueness in the interior of [0, oo)^ and complete the proof. □ 

B Annex : Classification Theorem for MBPRE 

We consider here a multitype branching process in random environment Z„ = {Zn^ : i = 
1, . . . , K) whose mean offspring matrix is denoted by ^ = A(w), where w = {wqjWi, . . .) is the 
environment. 

We introduce the extinction probability vector in environment w starting from one individual 
in habitat i: 

qi{w) = lim P(|Z„| =0\w, Z^^ = 1, Z^' = if i)- 

n— >-co 

Proposition B.l ( |20| . Theorems 9.6 and 9.10) Assuming that 

¥iq^{w) <l:i = l,...,K) = l or ¥{qiiw) = 1 : i = 1, . . . , K) = 1, (*) 

we have 

• //7 < 0, then the probability of extinction is equal to 1 for almost every w. 

• 7/7 = 1, then 

— either for every m>l, w-a.s., there exists 1 <i < K such that 
F{\Zra\ >l \ w, Zf = 1, = ^/j / i) = 0. 

— or qi{w) = 1 w-a.s.. 

Assuming that there exist integers N,L > such that P(V1 <i,j < K, {Ai\i ■ ■ ■ Ao)jj 7^ 0) = 1 
and |E(log(l — F{Zk = | Zq^^ = 1))| < 00, then (*) is satisfied and 

• 7/7 > 0, then w-a.s qi{w) < 1 for every i = I, • • • ,K, and 

P( lim log(|Z„|) = 7 I w^, 4'^ = 1, Zi = if j ^ i) = I - qi{w). 

^which states that two positive eigenvectors of a primitive matrix are colinear. It comes actually from the 
classical decomposition of A" using the maximum eigenvalue and the left and right eigenvector associated with. 
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Moreover thanks to |20^ Theorem 9.11], if ah the coefficients of the matrix A are positive and 
bounded, i.e., 

30 < c,c' < oo, c < inf Aij < sup Aij << c' a.s., 

^<i,j<K ' l<i,j<K 

then 

Zn = 0{\\An-i---Ao\\) a.s. 
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